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Abstract 

We compute instanton sizes and study correlation functions between instantons and 
monopoles in maximum abelian projection within SU{2) lattice QCD at finite tem- 
perature. We compare several definitions of the topological charge, different lattice 
actions and methods of reducing quantum fluctuations. The average instanton size 
turns out to be o" ~ 0.2 fm. The correlation length between monopoles and in- 
stantons is C ~ 0.25 fm and hardly affected by lattice artifacts as dislocations. We 
visualize several specific gauge field configurations and show directly that there is 
an enhanced probability for finding monopole loops in the vicinity of instantons. 
This feature is independent of the topological charge definition used. 



1 Introduction 



Classical gauge field configurations with non-trivial topology are believed to 
play an essential role in the confinement mechanism. In the scenario of the dual 
superconductor abelian monopoles condense leading to confinement. Large and 
interacting instantons could also produce confinement if they form an instan- 
ton liquid. At first sight these two pictures are distinct and the interesting 
question arises, whether instantons and monopoles are related to each other. 
Several groups investigated the relation between monopoles and instantons 
for semi-classical configurations [1,2]. We presented first evidence that those 
correlations also exist in realistic equilibrium configurations [3]. 

In this letter we compute instanton sizes from auto-correlation functions of 
topological charge densities. We calculate correlation functions between topo- 
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logical charge densities and monopole currents. For the extraction of the topo- 
logical charge we use field theoretical methods and the geometrical Liischer 
charge definition. In addition to the standard Wilson plaquette action we dis- 
cuss an eight parameter fix-point action [4]. Prom the auto-correlation func- 
tions of the topological charge instanton sizes can be estimated in smooth 
gauge-field configurations. Smoothing is done by cooling (for the Wilson ac- 
tion) and by means of blocking and reverse blocking (for the fix-point action) 
which is called constrained smoothing [5]. The correlation between monopoles 
and instantons is analyzed for the different topological charge definitions and 
after suppression of dislocations. Correlation functions are obtained both in 
the confinement and deconfinement phase of pure QCD. For direct insight into 
the local geometry of topological activity we visualize specific configurations 
by tools of computer graphics. 



2 Topological Chcirges and Monopoles 

There is no unique way to define a topological charge operator on the lat- 
tice. The field theoretical methods are straightforward discretizations of the 
continuum charge density: 



For our purposes we employ the plaquette and the hypercube method [6], 
consisting of a sum of products of link variables along two perpendicular pla- 
qucttcs or along a hypercube, respectively. The topological charge obtained 
has to be renormalized. A possibility to get rid of the renormalization con- 
stants is to use the method of cooling which reduces quantum fluctuations 
iteratively. Cooling, however, not only eliminates quantum fluctuations, but 
also small instantons and even large lattice instantons which die out if cooling 
is performed too intensely [7]. In particular we use the "Cabbibo-Marinari 
cooling method" with a cooling parameter oi 6 = 0.05. 

A way to overcome the problems associated with cooling is to use geometrical 
methods that interpolate the discrete set of link variables to the continuum in 
order to reconstruct the principal fibre bundle. We employ the locally gauge 
invariant Liischer charge definition for SU{2) [8]. A drawback of the geomet- 
rical methods is that for the Wilson action they are plagued by topological 
defects on the scale of the lattice spacing, dislocations. It it therefore necessary 
to smooth the gauge fields. To estimate the infiuence of dislocations on our 
results with the Wilson action, we use a simplified fix-point action that sup- 
presses dislocations [4]. To measure the sizes of instantons it is still necessary 
to suppress quantum fiuctuations. To this end constrained smoothing based on 
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renormalization group transformations was proposed as an alternative method 
to cooling which preserves long range physics and leaves instantons invariant 
[5]. 

In order to investigate monopole currents we project SU{2) onto its abelian 
degrees of freedom, such that an abelian U{1) theory remains [9]. This can be 
achieved by various gauge fixing procedures. We employ the so-called max- 
imum abelian gauge which is most favorable for our purposes. In our sim- 
ulations we subjected the configurations to 300 gauge fixing steps. For the 
definition of the monopole currents m{x, /i) we use the standard method [10]. 
After fixing the gauge the abelian parallel transporters u{x, //) are extracted 
and the color magnetic currents are computed: 

argM(n), (2) 

where u{^) denotes a product of abelian links u{x,iJ,) around a plaquette □ 
and f{x + fi, /i) is an elementary cube perpendicular to the /i direction with 
origin x + fi. From the monopole currents we define the local monopole density 
as p{x) = 4^E^|m(a;,/x)| . 



3 Correlation Functions 

Our simulations were performed on a 12^ x 4 lattice with periodic boundary 
conditions using the Metropolis algorithm. The obscrvables were studied for 
the Wilson action both in the confinement and the deconfinement phase at 
inverse gluon coupling /3 = A/g"^ = 2.25 (T/T^ = 0.88) and 2.4 (T/T^ = 1.29), 
respectively. For the fix-point action the inverse coupling ranged from (5 — 
4/c/2 = 1.50 (T/Tc = 0.83) to 1.8 (T/Tc = 1.73). For each run we made 100 
measurements, separated by 100 and 20 iterations for the Wilson action and 
for the fix-point action, respectively. 

The normalized auto-correlation functions (q'(O)g(r)) of the topological charge 
density are displayed in Fig. 1. In the case of the Wilson action they are pre- 
sented for the hypcrcubc definition (left) and for the Liischcr method (middle) 
for 0, 5, 20 cooling steps in the confinement phase at /3 = 2.25. Without cooling 
both auto-correlation functions are 5-peaked due to the dominance of quan- 
tum fiuctuations. They become broader with cooling refiecting the existence 
of extended instantons. The auto-correlation function of the hypercubc charge 
density is broader than that of the Liischer charge density, because the hyper- 
cube charge operator is more extended than Liischer's. The auto-correlation 
of the Liischer charge using a fix-point action is shown on the right-hand side 
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Fig. 1. Auto-corrclation functions of the topological charge density in the confine- 
ment phase using the hypercube definition (left) and the Liischer definition (middle) 
for 0, 5, 20 cooling steps for the Wilson action. Correlations for the Liischer charge 
definition before and after constrained smoothing for the fix-point action (right). 

of Fig. 1 in the confinement phase at /3 = 1.50 before and after constrained 
smoothing. It is again a (5-function for the original configurations and broad- 
ens after performing the smoothing procedure because quantum fluctuations 
are drastically reduced [11]. Here it is not necessary to finetune the cooling 
parameter and the number of cooling steps. 

To obtain an average instanton size we fitted the auto-correlation functions 
of the Liischer charge to the convolution f{x)= J Qa{t)Qa-{x — t) dt of the 
topological charge density Qa{x) = (^^i^)^ of a single instanton with 
size a. Such a fit is justified if instantons are dilute and well separated which 
is the case after 20 cooling sweeps or after constrained smoothing. The fitted 
instanton sizes a are displayed in physical units for various /S values in Table 1. 
The errors are in the range of 15%. For the Wilson action we used the 2-loop 
formula a{/3) = g{(3)/A with A = 6.6 MeV to obtain physical units. For the fix- 
point action we extracted the lattice spacing a at each /3 value from the string 
tension (440MeV)^ at T—0. Instantons show a tendency to become smaller on 
average with increasing temperature crossing the phase transition. However 
the results obtained with cooling seem less reliable due to the freedom in the 
cooling parameter. 

Table 1 

Instanton sizes a for different values of T/Tc computed from the auto-corrclation 
functions of the Liischer charge for the Wilson action after 20 cooling steps and for 
the fix-point action after constrained smoothing. 



T/T, 


0.83 


0.88 


0.92 


1.09 


1.29 


1.73 


a [fm] (Wilson action) 
a [fm] (Fix-point action) 


0.31 


0.20 


0.27 


0.21 


0.12 


0.10 
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Fig. 2. Correlation functions between the monopole density and the absolute value 
of the topological charge density for the hypercube definition (left) and the Liischer 
definition (middle) for 0, 5, 20 cooling steps for the Wilson action. Correlations 
using the Liischer definition and the fix-point action before and after constrained 
smoothing (right). The correlations are computed in the confinement phase (above) 
and in the deconfinement phase (below). 

As a measure for the local relation between abelian monopoles and instan- 
tons, we calculate the correlation functions {\q{0)\p{r)) between the absolute 
value of the topological charge density and the monopole density. They are 
displayed in Fig. 2 both in the confinement and the deconfinement phase, af- 
ter subtracting the cluster value and normalization. For the Wilson action the 
correlation functions are computed employing the hypercube method (left) 
and the Liischer method (middle) for several coohng steps. The shape of these 
correlations hardly changes under the influence of cooling and is essentially 
unaffected by the phase transition. Again the correlation functions with the 
hypercube charge are somewhat broader than those with the Liischer charge 
due to the different extensions of the operators. 

Liischer's charge together with the Wilson action is known to suffer from 
dislocations which might have a non-trivial correlation with monopoles [12]. 
To get rid of dislocations that may spoil the physical results, we use the fix- 
point action. The correlation functions of the Liischer charge distribution with 
the monopole density before and after constrained smoothing at T/T^ — 0.83 
and T/Tc = 1.09 are depicted on the right-hand side of Fig. 2. They turn 
out to be similar to those from the Wilson action and become slightly wider 
after smoothing. Therefore dislocations are not decisive for the non-trivial 
correlation found for the Wilson action. For all cases (Wilson action, fix-point 
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Table 2 

Ratio of spatial to time- like monopole density in the confinement and deconfinement 
for various cooling steps. 



Cool step 


P = 2.25 


p = 2.40 





0.996 


0.898 


5 


0.975 


0.571 


20 


0.268 


0.037 



action, cooling, constrained smoothing, confinement and deconfinement phase) 
the correlation lengths in lattice units are rather similar. The correlation length 
has a tendency to increase with temperature and was found to cover a range 
of C = 0.15 -0.35 fm. 

It has been reported that the ratio R = of space-like to time-like monopole 
densities decreases across the deconfinement phase transition and that it might 
serve as a reasonable order parameter [13]. We observe that the same quantity 
also decreases as a function of coohng which is displayed in Table 2. The 
drastic decrease yields some doubt on the quality of this quantity as an order 
parameter. For example after 20 cooling sweeps the string tension of SU (2) at 
(3 = 2.25 is still present, even though 55% of the monopole currents are static. 
However a.t (3 = 2.4 in the uncooled configurations only 27% of the monopoles 
are static but the string tension vanishes completely. 



4 Visualization 



We visualize the relation between instantons and monopoles by directly dis- 
playing clusters of topological charge and by drawing monopole loops in fixed 
time slices of specific configurations. For any value of the topological charge 
density q{x) > 0.01 a light dot and for q{x) < 0.01 a dark dot is plotted. The 
lines represent the monopole loops. 

In Fig. 3 we compare the results of the field theoretical methods and the 
Liischer method for the Wilson action on a single equilibrium gauge field after 
20 cooling sweeps. From left to right the plaquette, hypercube, and the Liischer 
charge distributions are plotted. The positions of the clusters of topological 
charge are the same for all methods. The points represent instantons or anti- 
instantons. In this particular configuration a monopole is found to wrap around 
the torus. 

Fig. 4 presents a cooling history of a time slice of a gluon field. The topolog- 
ical charge using the Liischer definition with the Wilson action is displayed 
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Plaquette definition Hypercube definition Liischer definition 




Fig. 3. Different definitions of the topological charge for a specific gluon configura- 
tion from Wilson action after 20 cooling sweeps. The instantons reside at the same 
places for all definitions and are surrounded by monopoles. In this configuration a 
monopole loop wraps around the torus. 

for cooling steps 0, 15, 20 and 25. Without cooling the topological charge dis- 
tribution cannot be identified with instantons due to quantum fluctuations. 
After 15-20 cooling steps one can assign instantons to clusters of topological 
charge. At cooling sweep 20 an instanton and an anti-instanton emerge. From 
cooling steps 35-40 they begin to approach each other and annihilate several 
steps later (not shown). Monopole loops also thin out with cooling, but they 
survive in the presence of instantons. In general, there is an enhanced proba- 
bility that monopole loops are present in the vicinity of instantons in all gauge 
field configurations. 



5 Conclusion 

We computed average sizes of instantons and correlation functions between 
instantons and monopoles for different topological charge definitions, for dif- 
ferent lattice actions, for cooled and constrained smoothed configurations. The 
auto-correlation functions between the topological charge density yield con- 
sistent results of instanton sizes for the Wilson and fix-point action. In the 
confinement phase the instanton size ranges from o = 0.2 — 0.3 fm whereas in 
the deconfinement phase a = 0.1 — 0.2 fm. The correlation functions between 
abelian monopoles and instantons are very similar for the geometrical and the 
field theoretical definition. They are hardly affected by cooling (Wilson action) 
or by constrained smoothing (fix-point action) and are qualitatively the same 
even across the deconfinement phase transition. The correlation length was 
found C = 0.15 — 0.25 fm in the confinement and ( = 0.25 — 0.35 fm in the 
deconfinement phase. 

We further calculated the local values of topological charges and monopole 
currents and directly displayed them with the help of computer graphics. Af- 
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Fig. 4. Cooling history for a specific gauge field configuration from Wilson action 
at a fixed time slice. The dots represent the topological charge distribution in the 
Liischer definition with \q{x)\ > 0.01. Positive (negative) charges are plotted with a 
light (dark) dot. The monopole loops correspond to lines. With cooling instantons 
evolve from noise accompanied by monopole loops in almost all cases. Note that 
time-like monopoles cannot be seen in these plots. 

ter a few cooling sweeps one observes clearly that instantons are accompanied 
by monopole loops. This correlation occurs on all (semi-classical) gauge field 
configurations considered. In a cooling history we demonstrated how instan- 
tons evolve from fluctuating gauge flelds and how they are surrounded by 
monopoles. Combining the above flnding that the correlations are rather in- 
sensitive under cooling or smoothing together with that of the 3D images, we 
conclude that the topological charge goes hand in hand with monopoles also 
in the original gauge fleld conflgurations. The demonstration of our simulation 
together with analytical investigations [1,14] might present a flrst indication 
of a deep relation between the topological structure of compact abelian and 
non-abelian gauge fleld theories. Since in abelian theories monopoles are re- 
sponsible for conflnement, and if such a relationship existed, this could be 
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accepted as a topological proof of quark confinement. 
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